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ABSTRACT 

The classical tests of general relativity — light deflection, time delay and perihelion shift 
— are applied, along with the geodetic precession test, to the five-dimensional extension of the 
theory known as Kaluza-Klein gravity, using an analogue of the four-dimensional Schwarzschild 
metric. The perihelion advance and geodetic precession calculations are generalized for the 
first time to situations in which the components of momentum and spin along the extra 
coordinate do not vanish. Existing data on light-bending around the Sun using long-baseline 
radio interferometry, ranging to Mars using the Viking lander, and the perihelion precession of 
Mercury all constrain a small parameter b associated with the extra part of the metric to be 
less than |6| < 0.07 in the solar system. An order-of-magnitude increase in sensitivity is possible 
from perihelion precession, if better limits on solar oblateness become available. Measurement 
of geodetic precession by the Gravity Probe B satellite will improve this significantly, probing 
values of b with an accuracy of one part in 10 4 or more. 

Subject headings: gravitation — relativity — solar system: general 
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Introduction 



There is now a substantial literature on the higher-dimensional extension of Einstein's general theory 



of relativity known as Kaluza-Klein gravity (Ovcrduin and Wesson 1997, Wesson 1999). There are several 
ways to test the theory, with perhaps the most straightforward involving the motion of test particles in 
the field of a static, spherically-symmetric mass like the Sun or the Earth. Birkhoff's theorem in the usual 



sense does not hold in higher dimensions (Bronnikov and Melnikov 1995, Schmidt 1997), so some question 
arises in identifying the appropriate metric to use for this problem. In the five-dimensional (5D) case 
(with one extra coordinate y = a; 4 ), most attention has focused on the soliton metric (Gross and Perry 



1983| , Sorkin 1983, Davidson and Owen 1985), which satisfies the 5D vacuum field equations, reduces to 
the standard four-dimensional (4D) Schwarzschild solution on hypersurfaces y = const, and contains no 
explicit y-dependence. The assumption of a vacuum in 5D is consistent with the spirit of Kaluza's idea, 
that 4D matter and gauge fields appear as a manifestation of pure geometry in the higher-dimensional 



world. The soliton metric has been generalized in various ways to incorporate time-dependence (Liu et al 



1993| ), y-dependence ( Billyard and Wesson 1996 ) and electric charge ( Liu and Wesson 1997| ), among other 
things (eg, Wesson and Liu 1998| ); see for review Overduin and Wesson (1997). We confine ourselves here 
to the original (two-parameter) soliton metric. 

The motion of test bodies in the gravitational field of the soliton can be studied using the familiar 
classical tests of general relativity (gravitational redshift, light deflection, perihelion advance and time 
delay), along with the geodetic precession test. Work done so far along these lines (Lim, Overduin and 
Wesson 1995; Kalligas, Wesson and Everitt 1995, hereafter "KWE") has demonstrated the existence of small 
but potentially measurable departures from the standard 4D Einstein predictions. In the present paper, we 
extend these earlier calculations in several ways, clarifying the physical meaning of the light deflection and 
time delay results for massless test particles and presenting new generalizations of the perihelion shift and 
geodetic precession formulas for massive ones. We take special care to compare our results to the latest 
experimental data in each case, obtaining new numerical constraints on the small parameter b associated 
with the extra part of the soliton metric. 



2. The Soliton Metric 



In what follows, lowercase Greek indices fi, v . . . will be taken to run over 0, 1, 2, 3 as usual, while 
capital Latin indices A,B,C... run over all five coordinates (0, 1, 2, 3, 4). Units are such that G = c = 1 
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except where stated otherwise. It is important to distinguish between the 4D line element (ds) and its 5D 
counterpart (dS), the two being related by 

dS 2 = ds 2 + g AA dy 2 . (1) 
To interpret an expression containing d/dS physically, one can always make the conversion 



d ds d / f dy\ 2 d 



,/s ,/S',/.s V 1 944 [dS J fa' (2) 



We emphasize in particular that d/dS ^ d/ds if dy/dS ^ 0. 



The soliton metric may be written (following Gross and Perry 1983, but switching to nonisotropic 
form, and defining a = 1/a, b = (3/ a and M = 2m) 

dS 2 = A a dt 2 - A- a - b dr 2 - A 1 - a ~ b x 

r 2 (d6 2 + sin 2 Ode/) 2 ) - A b dy 2 , (3) 

where A(r) = 1 — 2A//r, M is a parameter related to the mass of the object at the center of the geometry, 
and the constants a, b satisfy a consistency relation 

a 2 + ab + b 2 = 1 , (4) 

so that any two of M, a, b may be taken as independent metric parameters. We will treat b as the primary 
free parameter of the theory in what follows, noting that the 4D Schwarzschild metric is recovered (on 
hypersurfaces y = const) in the limit b —> (and a — > +1). In general, larger values of |6| will give rise to 
increasing departures from Einstein's theory, subject to the upper bound |6| < 2/\/3 ~ 1.15 imposed by 
equation (Q). Possible theoretical expectations for this parameter in the solar system and elsewhere are 
discussed further in sfi. 



3. Equation of Motion 

We proceed now with the analysis of experimental constraints. The Lagrangian for a test particle in 
the field described by the metric (||) is 

C = [A a i 2 - A- a - b r 2 - A 1 -"-* 1 x 

r 2 (6 2 + sin 2 0<P 2 ) - A b y 2 } 1/2 , (5) 
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where the overdot represents differentiation with respect to an affine parameter A along the geodesies. 
The Euler-Lagrange equations read 

d_ ( dC_\ _ dC_ _ 

We confine ourselves to orbits with 8 = ir/2 and 8 = 0, so that C becomes 

C = (AH 2 - A-°- b r 2 - A l - a - b r 2 tf - A b y 2 ) ^ . (7) 
We can identify three constants of motion 

C db 

where we have used the relation C = dS/dX. From these equations we find that 

—\ j_ Ar 2 ( — A 2 - 2a - b — A 2 - a - 2b 

d<p) +Ar ^ U 2 ~ ^ 2 

1 \2-a-b \ „4 



h2 A z - a - b j r 4 = . (9) 

The derivation here differs slightly from that of KWE, where £ = (dS/dX) 2 . Although the two approaches 
are physically equivalent, we have found that results are obtained more simply if the three constants of 
motion I, h, k are defined in terms of d/dS rather than d/dX (or d/ds). 



4. Light deflection 



Experimental upper limits on possible violations of local Lorentz invariance are extremely tight ( Will 



1991 ), so that we are justified in assuming that photons follow 4D null geodesies, ds — 0. The situation is 



not so clear with regard to the 5D line clement. However, it is economical to follow KWE and suppose that 
all particles follow ND null geodesies in N- dimensional gravity, whether massive or not.[] Proceeding on 



1 This assumption is supported by various lines of argument. In one version of 5D gravity, for example, the 
fifth coordinate y is related to rest mass m ( |Wcsson 198^ ), so that one has dS 2 = ds 2 + ga (G/c 2 ) 2 dm 2 . If all 
particles move on 5D null geodesies, then ds 2 = —g^(Gjc 2 ) 2 dm 2 . It then follows that ds = for photons, 
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this assumption, and substituting ds = dS = into equation (|l|), we get dy = also, so that £, h — > oo and 
k is undefined. The ratios Ijh and k/h are however well-behaved, and read 

i = r^A 2 ^- 1 $- = finite , 
h dcp 

k = r -2 A a+2b-l = Q _ (1Q) 

h dcp 

For self-consistency, therefore, the terms in k/h can be dropped from KWE equations (7), (8), (11) and (12). 
Equation (8) of that paper, in particular, reduces to 

and the definition of the parameter p, KWE equation (12), becomes just 

p=-(2-2a-b)L. (12) 

The photon's trajectory is deflected by an angle 

4M 

Sep = oj = h 2Mpr a , (13) 

r 

which agrees with KWE equation (18.1). At the closest approach to the central body, we have u = u a = l/r Q 
and du/dcj) = 0, so that equation ([ll] ) gives 

l jL=Al^ul = ±+0{e), (14) 

where e = M is a small parameter. Putting equations ( |l2| ) and ( |i"4| ) into equation (p^|), we find for the final 
light deflection angle 

(50= (4a + 26)— + 0{e 2 ) , (15) 
r a 

as in KWE equation (18.2), where however it is presented as a special case k = 0. We see here that 
equation ( [l5| ) is in fact entirely general for light deflection, and does not depend on any choice of k, which 
is in any case undefined when ds = dS = 0. 



which have m = const= 0. For massive particles with ds ^ 0, one expects variations in rest mass to, which 
are however below currently detectable levels, owing to the small size of the dimension-transposing constant 
G/c 2 (Overduin and Wesson 1997, 1998). Recent work on incorporating non-relativistic quantum theory into 
higher-dimensional gravity also strongly suggests that all test particles travel on ND null geodesies in the 



classical limit (3eahra 2000) 



-7- 



To obtain experimental constraints from the light deflection result, let us express equation (15) for the 
Sun in terms of the deviation A LD from the general relativity prediction <5</> GR , as follows 



W(l + A LD ) , 



(16) 



where (to first order in e): 



5<t> GIi ee 4M /r o , 

A LD ee a + 6/2-1 . (17) 

Using the consistency relation (^) we find 

a = -6/2±(l- 36 2 /4) 1 / 2 . (18) 

Theoretical and numerical work indicates that |6| <C 1 in the solar system (§|^), and our experimental limits 
bear this out. The negative roots of equation ( [l8| ) may also be ignored, as they are inconsistent with the 
limiting Schwarzschild case, and also imply the possibility of negative gravitational and/or inertial soliton 
mass ( Gross and Perry 1983| , Lim ct al. 1995, Overduin and Wesson 1997). We therefore take 

a= 1-6/2 - 36 2 /8 + 0(6 4 ) , (19) 

in the solar system, whereupon equation ( fl7| ) gives 

A LD = -36 2 /8 + 0(6 4 ) . (20) 

The best available constraints on A LD come from long-baseline radio interferometry, which implies that 
| A LD | < 0.0017 ( Robertson et al. 1991 , Lebach ct al. 1995 ). We therefore infer an upper limit 

|6|<0.07, (21) 

for the Sun. This could potentially be tightened by more than an order of magnitude using a proposed 
astrometric optical interferometer sensitive to departures from Einstein's theory of as little as |A LD | < 10~ 5 



(Rcascnbcrg and Shapiro 1986). 

It is important to bear in mind, however, that the parameter 6 characterizing the soliton metric (^) is 
not a universal constant of nature like G or c, but may in principle vary from soliton to soliton. Kaluza-Klein 
gravity as an alternative to 4D general relativity is therefore best constrained by the application of two or 
more tests to the same system. With this in mind we can use a recent measurement of light deflection by 



Jupiter, for which |A LD | < 0.17 (Ireuhaft and Lowe 1991), to obtain 

\b\ < 0.7, 



(22) 
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for that planet. It has also been proposed to measure light deflection by the Earth using the Hipparcos 



satellite, with an estimated precision of 12% (Gould 1993). Such a test would be sensitive to values of 
|6| < 0.6 for the Earth. The Gravity Probe B satellite should also be able to detect this effect by means of 



its guide star telescope, though with a somewhat lower precision (Adler 2000). 



5. Time delay 

The arguments in the previous section regarding the parameter k also apply to the time delay (or 
radar ranging) test, and circular photon orbits as well. That is, terms in k/h and k/t may be dropped from 
KWE equations (20-24) for radar ranging, and KWE equations (28-30) for circular orbits. The final results 
given there, however — equations (25) and (31) respectively — are correct. In fact, they hold not only for 
the special case k — 0, but quite generally. 

In particular, the excess round-trip time delay AT for signals emitted from Earth (at distance r e 
from the Sun) which graze the Sun (at nearest distance r ) and bounce off another planet (at r p ) may be 
calculated by setting k/£=0 in KWE equation (24) to obtain 



AT = AT GR (1 + A TD ) , (23) 



where (to first order in e) 



AT GR = 4M e 



In 



r£ - 7-2 



A TD = a + 6/2-1 

= -36 2 /8 + 0(6 4 ) . (24) 

We note that departures from 4D general relativity for time delay have exactly the same form as they do 
for light deflection. 

The best experimental constraint on time delay so far has come from the Viking lander on Mars, and 
gives | A TD | < 0.002 ( Rcascnbcrg et al. 1979| ). This leads immediately to the upper bound 



|6|< 0.07, (25) 

for the Sun, exactly the same as the limit obtained in the case of light deflection using long-baseline 
interferometry. 
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Keeping in mind that values of b can differ from soliton to soliton, however, it is possible that different 
physical setups could provide new information. For instance, one could attempt to measure b for the Earth 
by sending grazing signals from an orbiting satellite past our planet and bouncing them off the Moon; 



retroflectors left there by Apollo astronauts are routinely used for lunar laser ranging ( |Williams ct al. 1996|) . 
Substituting M e for M e and replacing r c ,r p and r with the appropriate distances, we find an expected 
excess time delay of order 400 ps using a satellite in geostationary orbit. This is well above the currently 



available resolution of ~ 50 ps (3amain et al. 1998). The feasibility of such a proposal would likely be 
limited by the weakness of the reflected signal. Better results might be obtained by active ranging between 
two orbiting satellites, or by statistical analysis of ranging data between two such satellites and an Earth 
station (the latter would however require excellent atmospheric modelling). 

In the same vein, one might attempt to measure b for the Moon by grazing it with signals from the 
Earth and bouncing them off the Viking lander on Mars. This might be done when Mars is at nearest 
approach (on the same side of the Sun as the Earth) to minimize signal contamination from the competing 
effect of the Sun. Substituting M m for M e in equation (p3|), however, and replacing r e ,r p and r with the 
appropriate distances, we find that the Moon's excess time delay (of order 10 ps) would be so short as to 
make this a daunting task at present. 



6. Perihelion advance 

We now switch our attention to massive test particles. In terms of a new variable u = 1/r, equation (^) 
becomes 

— \ + All 2 - f— 42-2a-6 _ jf_ A2-a-2b 

dd>) + Au U 2 ^ 2 

- ^ 2 ~ Q ~') = ■ (26) 

Differentiating with respect to <fi (and letting primes denote d/d<fi), we find that noncircular orbits (u' =/= 0) 
are governed by the following differential equation 

u" + (l + je)u = B + eB^u 2 + 0(e 2 ) , (27) 

where five new quantities have been introduced 



d= (2 - a - b) - l 2 {2 - 2a - b) 
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+k 2 (2 - a - 2b) , 
f = 2{2-a-b)(-l + a + b) 

+ 2£ 2 (-2 + 2a + 6)(-l + 2a + b) 

+ 2k 2 (2~a-2b)(-l + a + 2b) . (28) 

These expressions agree with KWE equations (32-36). (We have however chosen to relabel their e as /, for 
reasons that will become clear shortly.) 

The solution of the differential equation (f27j) is 



I = B + (l-2)CcoB{[l-e(l-2)],} 



u = - = B 
r 



C 2 f 7\ 2 
^-^ B ^2B( l -l) 

6 6B 



(l-|) 2 cos20 + O(e 2 ) , (29) 



where C is an integration constant. [This result differs slightly from KWE equation (37), where the factors 
of (1 — 7/2) 2 were omitted.] Equation ( p9[ ) can be written in a physically more transparent form by 
introducing two new quantities e and u) via 

(l-T)C = Be, l- e (l-l) sw . (30) 

With these definitions, we find that 

u = - = B(l + ecoswcf)) + -eB 2 \-e 2 cos2uj<p 
r 2 

Xl +0(e 2 ), (31) 



6(l-7+2 P/ 



where e = M is a small parameter as before, and 



uj = 1 - 3eB ( 1 + ) + 0(e 2 ) . (32) 



6d, 

The first term on the right-hand side of equation ( |3l| ) is of order e°, and shows explicitly the elliptical shape 
of the orbit. This is then modified by the second term, of order e . Note that e is just the eccentricity of 
the ellipse. The angular shift between two successive perihelia is given by 

QnM 2 d ( f 



in agreement with the final result (38.1) of KWE. It should be emphasized that the angular momentum h 
is not in general the same quantity in 5D as it is in 4D. In particular, putting equations (||) and (||) into 
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the second of equations (H), we find 



A 



\-a-bJl ' 



dS 



A l-a-b r 2^_^ 1 + A - bk 2 

ds 



/i (4D) v / l + A- b k 2 



(34) 



If k 0, therefore, it follows that h ^ h 



(4D)- 



To eliminate h from equation (|33|), let us consider the points along the orbit where r takes its 
minimum value r_ and maximum value r + respectively. From inspection of equation ( |3l| ) we see that 
r_ = B~ l (\ + e) _1 + 0(e) at uxp = and r + — B~ l (\ — e) _1 + 0(e) at uj(f> = tt. The semimajor axis a Q of 
the ellipse is then 

0(e), (35) 



so that 



D 



1 



/i 2 Oo(l - e 2 ) 



0(e) 



(36) 



h 2 =e(l-e 2 )a o d + 0(£ 2 ) 



Substituting equation ( p6[ ) into equation (33), we find 

50 



6ttM ( ! \ 



(37) 



(38) 



Only one term in this result remains physically obscure, and that is the ratio f/d. This is given in terms of 
I and k by the definitions (Eq). The latter two constants are related by equation (E6I) as follows 



t 2 = h 2 



$M +Au 2 



A 



2a+b-2 , ^ 2 y 4 a - b 



(39) 



Since h 2 is of order e 1 by equation (^), while u and it' are of order e° by equation (^Tj), it follows from 
equation ( p)| ) that £ 2 = 1 + fc 2 + 0(e). Using the definitions (pj|), we therefore obtain 

/ 26 /c 2 6(a-6)/3 rt/ . 

6a 3 a + k z (a — 6) 



so that the final perihelion precession angle (38) becomes 



6nM 



2 k 2 (a- 6)6/3 
fl+ 3 a + P(a-6) 



0(e 2 



(41) 



r ao(l-e 2 ) 

This represents the generalization of KWE equation (38.2) to cases in which k ^ (and eccentricity 
e / 0). In the special case 6 = (and a = +1), for which the metric (|J) reduces to Schwarzschild form on 
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hypersurfaces y = const, it is interesting to note that one recovers the standard 4D general relativity result, 
regardless of the value of k. In this limit, therefore, the perihelion shift test is insensitive to the momentum 
of the test body along the extra coordinate. And in general, one must choose a soliton with 6 ^ in order 
to distinguish experimentally between test particles with different values of k. 

As usual, let us parametrize our result in terms of the departure from 4D general relativity so that 

5<t> = ty QR (l + A PP ) , (42) 



where (to first order in e): 



QttM 



Ml - e 2 ) ' 



2, fc 2 (a- 6)6/3 

A PP = a+-b+ 2 l2( h\ - 1 ' 43 
3 a + fir (a — b) 

Theoretical work indicates that k, which is a measure of momentum along the fifth dimension, is related to 



the charge-to-mass ratio of the test body ( Wesson and Liu 1997 ). For a planet such as Mercury, we may 



take k = 0. Putting equation ( |19[ ) into equation (|43|), we therefore have 

A PP = 6/6 - 36 2 /8 + <3(6 4 ) . (44) 

Perihelion precession is thus a potentially more sensitive probe of higher-dimensional gravity than either 
light deflection or time delay, in that it depends on the first, as well second order in b. 

Unfortunately, however, this increased sensitivity is offset in the case of Mercury's orbit about the Sun 
by uncertainty in the solar oblateness. The latter introduces a new term £J 2 inside the brackets on the 



right-hand side of equation (42), where £ = i?g/2M a o (l — e 2 ) and J 2 is the solar quadrupole moment 
( Campbell et al. 1983 ). Dividing through by the orbital period T, we may therefore write for the rate of 
perihelion advance (to order 6 3 ) 

Aw = y- Aw GR (l + £ J 2 + 6/6 - 36 2 /8) , (45) 

where Aw GR = S(p GR /T — 42.98 arcsec/century. The observed value of Mercury's perihelion precession 
rate is quite close to this value, Auj = 43.11 ± 0.21 arcsec per century (Shapiro, Counselman and King 
1976). Experimental data on J 2 is a good deal more controversial and has ranged over two orders of 



magnitude, from a maximum value of (23.7 ± 2.3) x 10 6 (Dickc and Goldenberg 1967) to a minimum of 
(0.17 ± 0.02) x 10~ 6 ( puvall et al. 1984| ). One straightforward least -sqares fit to a number of published 



measurements leads to intermediate value of J 2 = 5.0 x 10 6 , which however implies a general relativistic 
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precession rate more than two standard deviations away from that observed (Campbell ct al. 1983). Such 
a discrepancy could be explained in the context of higher-dimensional gravity by modelling the Sun as a 
soliton with b = —0.062. This is just consistent with the constraint \b\ < 0.07 from light deflection (§[|) and 
time delay (§|J), which is intriguing since these tests probe somewhat independent aspects of relativistic 
gravity. Improved experimental data relating to any of the three tests would be of great interest. 

Conservative limits on b from perihelion precession may be obtained by quoting the results of a recent 
review in which all available data (to 1997) have been combined to give a weighted mean value for the solar 



oblateness of J2 = (3.64 ± 2.84) x 10 6 ( Rozclot and Rosch 1997 ). Using this uncertainty range, together 
with that in the observed value of Aw for Mercury's orbit, we hnd that 

b = -0.03 ±0.07, (46) 

for the Sun. This is consistent with the bounds obtained from light deflection and time delay. Sensitivity of 
the perihelion precession test to the value of b could be improved by an order of magnitude if better data on 
J2 were to become available; the proposed ASTROD mission, for example, might measure this parameter 



to an accuracy of 5 x 10 8 (Ni 199£ 



7. Geodetic effect 

We now move on to consider spinning massive test particles with velocity 5-vectors u = dx j dS and 
spin 5-vectors S . The motion of these objects is governed by three central equations; namely, the geodesic 
equation 

f^ + f c AB u A u B =0, (47) 
^+f c AB S A u B = 0, (48) 



the parallel transport equation 



and the orthogonality condition 

u c S c = 0. (49) 

Here T AB refers to the 5D Christoffel symbol for the metric (Q). This is defined in exactly the same 
manner as the usual 4D Christoffel symbol, with indices running over five values instead of four (see KWE, 
Appendix Al for details^). 



2 There are some minor typographical errors in this appendix, which we note briefly here. The factors 
of (1 - 2M)/r in equations (A2.2), (A2.6) and (A2.7) should read 1 - 2M/r. The same thing applies to 
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In order to simplify the problem, we follow KWE in assuming that the test particle moves in a circular 
orbit with 8 = tt/2, r = r D and 9 = r = 0. Its velocity u c may then be expressed as follows in terms of the 
constants of motion £, h and k, as given by equations (|^) 

dx c 



u 



c = (lA~ a , 0, 0, hr~ 2 A a+b - 1 , kA- b ) . (50) 



dS 

From the metric (|^), we have 

1 = A a (u a ) 2 - A^^rl (u 3 ) 2 - A" (u 4 ) 2 , (51) 
which, with equation (p0[), implies 

£ 2 - h 2 r- 2 A 2a+b - 1 - k 2 A a - b — A a = . (52) 
It may be shown that the motion of the test body as given by equations (pO) and (p2|) is geodesic in the 



sense of equation (47) 



We now propose to generalize the treatment of KWE by leaving the extra component S 4 of spin 
unrestricted, rather than setting it to zero. In fact, writing explicitly S c = (S° , S 1 , S 2 , S 3 , S 4 ), we find 
that the orthogonality condition ( p9| ) imposes the following restriction on the spin components 

£S° - hS 3 - kS 4 = , (53) 

so that S will not vanish in general, if the parameter k is well-defined. 

We now proceed to solve the parallel transport equation dig) , taking one value of the index C at a 
time. To begin with, the C = 2 component gives 

S 2 = — = const , H 2 = const . (54) 
r a 

(Note that, due to our choice of coordinates, S ^ 1 and S 4 are dimensionless while S 2 and S 3 have units 
of inverse length.) Defining a new function g = g(S) of the 5D proper time, we may write without loss of 
generality 

S° = H g, H = const . (55) 



equations (57) and (58) in the main body of KWE. Also, the exponents —(1/2) and 1/2 in equations (57) and 
(58) should be switched, in agreement with equations (A2.7) and (A2.2) respectively. These discrepancies 
do not affect any of the other equations or conclusions reported in KWE, and do not appear in the new 
reference book on Kaluza-Klein gravity by Wesson (1999). 
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The C = component of equation (|48|) then reads 
The C = 4 component, meanwhile, takes the form 



H r 2 Aa+1 

aiM ° dS 



(56) 



S 4 = H 4 g + K 4 , 

H 4 = ^H A a - b 
at 

Ka = const , 



= const , 



(57) 



where we have used equation (|56|) . In a similar way, the C = 3 component of equation ( |48| ) gives 

S 3 = ^g + K 3 , 

K 3 



>o 

hH 
aiM 
const . 



1- (l + a + b — 



= const , 



(58) 



We now solve the C = 1 component of equation (Mq), assuming for simplicity that K3 = K 4 = 0. Using 



equations (J50J) , (|55|), ( |57| ) and (J5SJ) , we find 



rfS 1 
~d~S~ 



_E2E A -+b-l! a 2 e 2_ b 2 k 2 A a- b 

atri I 



A/ 2 



1- (l + a + b) 



M 



A 



2a+fe-l I 



Differentiating equation (pq) with respect to S, meanwhile, gives 



dS 1 
~dS~ 

Equating these two expressions, we obtain 



J^_ r 2 A a+i f_9_ 
aiM dS 2 



dS 2 



= -n 2 g 



where 



h 2 

r 4 



l-(l + a + 6) — 



A 



2a+b-l 



(59) 



(60) 



(61) 



(62) 
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The general solution of equation ( |61| ) is g(S) = K\ sin(f75 f ) + K2 cos(f25). We choose K\ — 1 and K2 = 
for simplicity. The spin components are then given by 



S = H frin(fiS) , S 1 = #1 cos(flS) , 
S 4 = H 4 sm(nS) , 



(63) 



where Hi and Hi are arbitrary constants and 

alM 



Hn = - 



m = — —A 



Ha = — 



h 

6fcM 



+~ Hi 



ia+fc-2 



1- (1 + + 6) 



M 



Hi 



r 2 il 



A- b - l Hi . 



(64) 



The spatial part of S c is thus seen to rotate in the plane of the orbit with angular speed Q. Substituting 
these results into equation (m) yields 



A 



2a+b-l 



Mr D 
bk 2 A a - b = 



1- (I + + 6) — 



Solving simultaneously with equation (p3), we obtain for the constants of motion 



(65) 



2 A-b 



A a <l + k*A 



l-a-b 



M 



a+(a- b)k 2 A 



2 A-b 



1 - (1 + 2a + b)M/r 



h 2 = Mr n A 



2 A-b 



a + (a - b) k 2 A 



1 - (l + 2a + 6)M/r _ 
These expressions can be written in terms of a small parameter e = M as usual 

bk 2 \ M~ 



(66) 



I 2 = (1 + k 2 ) 



1- a 



1 + fc 2 



+ 0(e 2 ), 



h 2 = Mr„ [a + (a - 6) fc 2 ] {1 + [(4a + 36-1) 



+ 



26(a - 6) k 2 



AI 



+ 0(e 3 ) 



a + (a — 6) fc 2 

With the aid of equation (p2), we then find for the angular speed of the spin vector 



(67) 



b(a~b)k 2 1 



a + (a - b) k 2 



0{e 2 



(68) 
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This quantity is not the same as the test body's orbital angular speed, which is given in terms of the 5D 
proper time dS as 



W S f s = hr- 2 A a+b - 1 



(69) 



A (a+b-i)/2 / a+(a-b)k 2 A 



2 A-b 



1 - (l + 2a + b)M/r ' 
where we have used equations (|3|) and (|66|). In terms of e 

M / 3 — 6 



[a + (a - 6) /c 2 ] M 



6(a - 6) k 2 



1 



(70) 



a + (a — 6) A: 2 / 

It is precisely the excess of f2 over ui that gives rise to the geodetic effect. 

Suppose the spin vector S c is initially oriented in the radial direction; ie, if 2 = at S — 0. During one 
orbit, the test body's angular displacement <f> goes from to 2ir, so that SS = 2-k/uj. In the same period, 
S 3 goes from its initial value of zero at S — to its final value at proper time S. To first order in e, the 
spin vector has advanced through an angle 

^_r o [S 3 (S)~S 3 (0)] 



S\0) 
Hi V w 



+ 0(e 2 ), 
0(e 2 ), 



-2n 



n 



-1 ) 



(71) 



where we have used equations (p4[), (|6J) and (pq). Combining equations ( pq ) and (|7^), we find that 

+ 0(e 2 ), (72) 



oj 2r Q 



2 fc 2 6(a-6)/3 
fl+ 3 a+(o-6) fc 2 



so that the geodetic precession angle can finally be expressed as follows in terms of its deviation from the 
prediction of 4D general relativity 

5(b = 5(b GR (1 + A GP ) , (73) 



where (to first order in e) 



3irM/r 
2 

a+-b, 



k 2 {a-b)b/3 
a+k 2 (a - b) 



(74) 
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This represents the generalization of KWE equation (66) to cases in which S 4 ^ 0. Deviations from 4D 
general relativity have exactly the same form for geodetic precession as they do for perihelion precession. 
Taking k = and using equation (|l^), as in ^6[ we find that 

A GP = b/6-3b 2 /8 + 0(b 4 ) . (75) 

Like the perihelion shift, geodetic precession depends on b to first as well as second order, and is thus a 
potentially more sensitive probe of the theory than either light deflection or time delay. 

The Gravity Probe B satellite, currently scheduled for launch in early 2001, has been designed to 
measure deviations from 4D general relativity with a precision of better than |A GP | < 2.5 x 10~ 5 (Buchman 



et al. 1996 ). Using equation (|75|), we find that this corresponds to a sensitivity to values as small as 

\b\ < 1 x 10~ 4 , (76) 

or better for the Earth — a constraint some five hundred times stronger than any other solar system bound 
obtained to date, and five thousand times stronger than the only other Earth-based test (light deflection 
using Hipparcos; §||). 

We conclude this section by noting that a complementary analysis of geodetic precession has been 
carried out for a static, spherically-symmetric 5D metric different from that given by equation ([|), one 
in which the fifth dimension is flat (Mashhoon, Liu and Wesson 1994; Mashhoon, Wesson and Liu 1998). 
The inclusion of spin is of special importance in this case since the classical tests (based on the equations 
of motion) alone cannot discriminate between 4D and 5D effects. The geodetic precession rate has been 



computed, and differs from the 4D Einstein value in the weak-field, low velocity limit (Liu and Wesson 



1996). A preliminary interpretation of the discrepancy indicates, however, that it is likely to be somewhat 



below the threshold of detection by Gravity Probe B (Overduin and Wesson 1998) 



8. Discussion 

Having obtained upper limits on |6| of order 0.07 (and possibly 10~ 4 ) from experiment, we consider 
here the range of values that might be expected for this parameter on theoretical grounds. These turn out 
to be small (perhaps of order 10 -8 to 10 -2 ) in the solar system, but could be larger (of order 0.1) in larger 
systems such as clusters of galaxies. 

These estimates are based on the fact that the soliton's effective 4D mass is not concentrated at a point, 
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like that of a black hole, but has instead a finite (though sharply peaked) density profile whose steepness 
depends on the metric parameters (Liu and Wesson 1992, Wesson and Ponce de Leon 1994| ). Quoting the 
latter authors, but replacing their metric parameters a, e, k (due to Davidson and Owen 1985) with our 
M, a, b via M = 2/3, a = ek and b = — e, we find for the density of the soliton 

n , , -abM 2 /^ (l~M/2r\ 2(a+b) 



[1 - (M/2r) 2 ] 4 \l + M/2r 
Pressure is given by p = p/3, so that the matter described by equation (fr^) could be radiationlike, or 
composed of ultrarelativistic particles such as neutrinos. Total gravitational mass (as deduced from the 
asymptotic form of the metric) is M g — aM, so it is clear that b must be negative for positive density. 
Numerical analysis further reveals that the mass of the soliton is increasingly concentrated at small r as 
|6| approaches zero, and that the 4D Schwarzschild limit (b — 0) can in fact be viewed as a maximally 



compressed soliton (Wesson and Ponce de Leon 1994). Physically, this means that solar system bodies, 
which (viewed as solitons) are essentially point masses, are likely to be associated with very small values of 
\b\. 



To attach some numbers to these qualitative remarks, we make use of equation ([19|) and consider the 
weak-field (r 3> M/2), small-6 limit, in which 

p{r)^~bGM 2 g /(^c 2 r A ) , (78) 

where we have reverted to physical units. Equation (|T8|) allows us to associate ranges of 6-values with 
solitons of mass M g , if the density p can be estimated at some radius r. It has, for instance, been suggested 



(eg, Freese 1986 , Gould 1992 ) that relativistic hot dark matter in the form of massive neutrinos could be 
trapped inside the Earth. Krauss et al. (1986) have derived one possible density profile for such particles, 
assuming that equilibrium is established between those undergoing capture, annhilation, and escape from 
the Earth's gravitational potential. We do not attempt to fit our equation ( f78| ) to this profile at all radii, 
but merely take the predicted neutrino density at the Earth's surface as illustrative. From Fig. 2 of Krauss 
et al. (1986), the expected escape rate for 10 GeV neutrinos is 2 x 10 16 s _1 , which translates into a density 
at the Earth's surface of p(R®) = 3 x 10 -20 kg m~ 3 (about 50 times the canonical local halo dark matter 
density of 5 x 10~ 22 kg m~ 3 ). If we suppose that this is rather associated with solitonic matter making 
up some fraction £ of the Earth's total mass (M g = (Mq), then equation (|7Sl) gives b = —4 x 10~ 14 C~ 2 . 
For dark matter of this kind to be significant, b must be small for solar system bodies. With £ ~ 10~ 3 , 
for example, we have b ~ —4 x 10~ 8 , while £ ~ 10 -6 would correspond to b ~ —0.04. These numbers 
are consistent with the experimental limits obtained in §§ |J - [t] above. It may be possible to constrain 
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the theory more tightly by looking at violations of the weak equivalence principle by solar system bodies 
( pverduin 2000| ). 

On larger scales, systems such as galaxies and clusters of galaxies are suspected by many to harbor 
significant amounts of relativistic hot dark matter. We take here as an example a recent numerical 



simulation (Kofman et al. 1996) in which light (2.3 eV) neutrinos make up 20% (by mass) of a cluster 
whose total mass M T = 6 x 1O 14 M0. Fig. 3 of this paper shows a typical neutrino density of p(r) w 200p c 
at r — 0.03 Mpc, where p c = 2 x 10~ 26 ft.Q kg m~ 3 is the critical density. If this were instead attributed to 
solitonic dark matter of total mass M g = QM T , then the latter would have b = — 0.01£~ 2 by equation (|7q) , 
where we have taken ho — 0.65. If all the hot dark matter were solitonic (£ = 0.2), then |6| could be as 
large as 0.3. These values are illustrative only, since density profiles of hot dark matter in clusters are likely 
somewhat shallower than that indicated by equation (J7^)Q Nevertheless they establish that values of |6| in 
galaxy clusters might in principle be significantly larger than those in the solar system, and this encourages 
us to speculate that stronger tests of higher-dimensional gravity might be carried out using the excellent 
observational data now available on gravitational lensing by these objects. 



9. Conclusions 

We have re-examined the classical tests of general relativity, as well as the geodetic precession test, 
when Einstein's theory is extended from four to five dimensions. The physical meaning of previous 
calculations for light deflection and time delay have been clarified physically, and the restriction of zero 
momentum and/or spin along the extra coordinate that characterized the earlier calculations of perihelion 
shift and geodetic precession has been lifted. 

Our results show that Kaluza-Klein gravity remains consistent with experiment. The free parameter 
of the theory, however, is increasingly constrained to small values. Thus, data on light deflection, radar 
ranging to Mars and the perihelion precession of Mercury all imply a value of \b\ < 0.07 for the Sun. 
Improved data on solar oblateness should improve the sensitivity of the perihelion precession bound by as 
much as an order of magnitude. And the upcoming launch of Gravity Probe B will allow us to measure 
values of |6| for the Earth with an accuracy of one part in 10 4 or better. 



1 Density profiles with p oc r 4 at large r have however been discussed in other contexts, such as elliptical 



galaxies (Jaffc 1983, de Zeeuw 1985, Hernquist 1990) 



- 21 - 



We thank P. S. Wesson for comments. H. L. acknowledges the support of the National Natural Science 
Foundation of China (grant no. 19975007). J. O. thanks R. I. Bush and W.-T. Ni for comments on solar 
oblateness and radar ranging, and acknowledges the support of the National Science and Engineering 
Research Council of Canada. He also expresses his gratitude for the hospitality of C. W. F. Everitt and the 
theory group at Gravity Probe B, Stanford University, where part of this work was carried out. 



- 22 - 

REFERENCES 

Adlcr, R. J. 2000, Gravity Probe B preprint, Stanford University 
Billyard, A. and Wesson, P. S. 1996, Phys. Rev., D55, 731 
Bronnikov, K. A. and Melnikov, V. N. 1995, Gen. Rel. Grav., 27, 465 

Buchman, S. et al. 1996, in proc. 7th Marcel Grossmann Meeting on General Relativity, ed. R. T. Jantzen 
and G. M. Keiser (Singapore, World Scientific), p. 1533 

Campbell, L., McDow, J. C, Moffat, J. W. and Vincent, D. 1983, Nature, 305, 508 

Davidson, A. and Owen, D. A. 1985, Phys. Lett., B155, 247 

de Zeeuw, T. 1985, Mon. Not. R. Astron. Soc, 216, 273 

Dicke, R. H. and Goldenbcrg, H. M. 1967, Phys. Rev. Lett., 18, 313 

Duvall, T. L., Dzicmbowski, W. A., Goodc, P. R., Gough, D. O., Harvey, J. W. and Leibacher, J. W. 1984, 
Nature, 310, 22 

Freese, K. 1986, Phys. Lett., 167B, 295 

Gould, A. 1992, ApJ, 388, 338 

Gould, A. 1993, ApJ, 414, L37 

Gross, D. J. and Perry, M. J. 1983, Nucl. Phys., B226, 29 

Hernquist, L. 1990, ApJ, 356, 359 

Jaffe, W. 1983, Mon. Not. R. Astron. Soc, 202, 995 

Kalligas, D., Wesson, P. S. and Evcritt, C. W. F. 1995, ApJ, 439, 548 

Kofman, L., Klypin, A., Pogosyan, D. and Henry, J. P. 1996, ApJ, 470, 102 

Kolb, E. W. and Turner, M. S. 1990, The Early Universe (Reading, Addison- Wesley) , p. 441 

Krauss, L. M., Srednicki, M. and Wilczek, F. 1986, Phys. Rev., D33, 2079 



- 23 - 

Lebach, D. E. ct al. 1995, Phys. Rev. Lett., 75, 1439 

Lim, P. H., Ovcrduin, J. M. and Wesson, P. S. 1995, J. Math. Phys., 36, 6907 

Liu, H. and Wesson, P. S. 1992, J. Math. Phys., 33, 3888 

Liu, H. and Wesson, P. S. 1996, Class. Quant. Grav., 13, 2311 

Liu, H. and Wesson, P. S. 1997, Class. Quant. Grav., 14, 1651 

Liu, H., Wesson, P. S. and Ponce dc Leon, J. 1992, J. Math. Phys., 34, 4070 

Mashhoon, B., Liu, H. and Wesson, P. S. 1994, Phys. Lett., B331, 305 

Mashhoon, B., Wesson, P. and Liu, H. 1998, Gen. Rel. Grav., 30, 555 

Ni, W.-T. 1998, J. Jpn. Soc. Microgravity Appl., Vol. 15, Suppl. II, 66 

Ovcrduin, J. M. 2000, Phys. Rev. D., submitted 

Overduin, J. M. and Wesson, P. S. 1997, Phys. Rep., 283, 303 

Overduin, J. M. and Wesson, P. S. 1998, in Current Topics in Mathematical Cosmology, ed. M. Rainer and 
H.-J. Schmidt (Singapore, World Scientific), p. 293 

Reasenberg, R. D. and Shapiro, I. I. 1986, in Relativity in Celestial Mechanics and Astrometry, ed. J. 
Kovalevsky and V. A. Brumberg (Dordrecht, Reidel), p. 383 

Reasenberg, R. D. et al. 1979, ApJ, 234, L219 

Robertson, D. S., Carter, W. E. and Dillingcr, W. H. 1991, Nature, 349, 768 
Rozelot, J. P. and Rosch, J. 1997, Solar Phys., 172, 11 
Samain, E. et al. 1998, Astron. Astrophys. Suppl. Ser., 130, 235 
Schmidt, H.-J. 1997, Grav. Cosm., 3, 185 
Seahra, S. S. 2000, Phys. Rev. D., submitted 

Shapiro, I. I., Counselman, C. C. and King, R. W. 1976, Phys. Rev. Lett., 36, 555 



-24- 

Sorkin, R. D. 1983, Phys. Rev. Lett., 51, 87 
Trcuhaft, R. N. and Lowe, S. T. 1991, AJ, 102, 1879 
Wesson, P. S. 1984, Gen. Rcl. Grav., 16, 193 

Wesson, P. S. 1999, Space, Time and Matter (Singapore, World Scientific) 

Wesson, P. S. and Liu, H. 1997, Int. J. Thcor. Phys., 36, 1865 

Wesson, P. S. and Liu, H. 1998, Phys. Lett., B432, 266 

Wesson, P. S. and Ponce de Leon, J. 1994, Class. Quant. Grav., 11, 1341 

Will, C. M. 1993, Theory and experiment in gravitational physics (Cambridge University Press, Cambridge), 
§2 

Williams, J. G., Newhall, X. X. and Dickey, J. O. 1996, Phys. Rev., D53, 6730 



This manuscript was prepared with the AAS L^TjtX macros v4.0. 



